EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 



CERN-PH-EP /2006-040 
20 December 2006 



Spin asymmetry Af 

and the spin-dependent structure function gf 
of the deuteron at low values of x and Q 2 



The COMPASS Collaboration 



Abstract 

We present a precise measurement of the deuteron longitudinal spin asymmetry Af 
and of the deuteron spin-dependent structure function gf at Q 2 < 1 (GeV/c) 2 and 
4 • 10~ 5 < x < 2.5 • 1CT 2 based on the data collected by the COMPASS experiment 
at CERN during the years 2002 and 2003. The statistical precision is tenfold better 
than that of the previous measurement in this region. The measured Af and gf are 
found to be consistent with zero in the whole range of x. 

Keywords: inelastic muon scattering; spin; structure function; Al; gl; low x; low 
Q2. 



(Submitted to Physics Letters B) 



The COMPASS Collaboration 



E.S. Ageev 24) , V.Yu. Alexakhin 8) , Yu. Alexandrov 18 ) , G.D. Alexeev 8) , A. Amoroso 29 ), 
B. Badelek 30 ), F. Balestra 29 ), J. Ball 25 ), G. Baum 1 ), Y. Bedfer 25 ), P. Berglund 13 ), 
C. Bernet 25) , R. Bertini 29 ), R. Birsa 28) , J. Bisplinghoff 3) , P. Bordalo 15 - a) , 
F. Bradamante 28) , A. Bravar 16 ), A. Bressan 28 ), E. Burtin 25 ), M.P. Bussa 29) , 
V.N. Bytchkov 8 ), L. Cerini 28 ), A. Chapiro 27 ), A. Cicuttin 27 ), M. Colantoni 29 ' b ), 
A. A. Colavita 27 ), S. Costa 29 ), M.L. Crespo 27 ), N. d'Hose 25 ), S. Dalla Torre 28 ), 
S.S. Dasgupta 6 ), R. De Masi 20 ), N. Dedek 19 ), O.Yu. Denisov 29 ' c ) , L. Dhara 7 ), V. Diaz 
Kavka 27 ), A.M. Dinkelbach 20 ) , A.V. Dolgopolov 24 ), S.V. Donskov 24 ), V.A. Dorofeev 24 ), 
N. Doshita 21 ), V. Duic 28 ), W. Diinnweber 19 ) , J. Ehlers 12 ' 16 ), PD. Eversheim 3 ) , 
W. Eyrich 9 ), M. Fabro 28) , M. Faessler 19 ), V. Falaleev 11 ), P. Fauland 1 ), A. Ferrero 29) , 
L. Ferrero 29 ), M. Finger 22 ), M. Finger jr. 8 ), H. Fischer 10 ), J. Franz 10 ), J.M. Friedrich 20 ) , 

V. Frolov 29 ' c ), U. Fuchs 11 ), R. Garfa gnini 29 ), F. Gautheron 1 ), O.P. Gavrichtchouk 8 ) , 
S. Gerassimov 18 ' 20 ), R. Geyer 19 ), M. Giorgi 28 ), B. Gobbo 28 ), S. Goertz 2 ' 4 ), A.M. Gorin 24 ), 
O.A. Grajek 30) , A. Grasso 29) , B. Grube 20) , A. Griinemaier 10 ) , J. Hannappel 4) , D. von 
Harrach 16 ), T. Hasegawa 17) , S. Hedicke 10) , F.H. Heinsius 10) , R. Hermann 16 ), C. Hefi 2) , 

F. Hinterberger 3 ) , M. von Hodenberg 10 ) , N. Horikawa 21,d ), S. Horikawa 21 ), 
R.B. Ijaduola 27 ), C. Ilgner 19 ), A.I. Ioukaev 8 ), S. Ishimoto 21 ) , O. Ivanov 8 ), T. Iwata 21 ' 6 ), 
R. Jahn 3) , A. Janata 8 ), R. Joosten 3 ), N.I. Jouravlev 8 ), E. KabuB 16) , V. Kalinnikov 28 ) , 
D. Kang 10 ), F. Karstens 10 ), W. Kastaun 10 ), B. Ketzer 20 ), G.V. Khaustov 24 ), 
Yu.A. Khokhlov 24 ), N.V. Khomutov 8 ), Yu. Kisselev 1 - 2 ), F. Klein 4 ), S. Koblitz 16 ), 
J.H. Koivuniemi 13 ) , V.N. Kolosov 24 ), E.V. Komissarov 8 ) , K. Kondo 21 ), 
K. Konigsmann 10 ) , A.K. Konoplyannikov 24 ) , I. Konorov 18 ' 20 ) , V.F. Konstantinov 24 ) , 
A.S. Korentchenko 8 ) , A. Korzenev 16,c ) , A.M. Kotzinian 8 ' 29 ), N.A. Koutchinski 8 ) , 
K. Kowalik 30 ), N.P Kravchuk 8 ), G.V. Krivokhizhin 8 ) , Z.V. Kroumchtein 8 ) , R. Kuhn 20 ), 
F. Kunne 25 ), K. Kurek 30 ), M.E. Ladygin 24 ), M. Lamanna 11 ' 28 ), J.M. Le Goff 25 ), 
M. Leberig 11 - 16 ), J. Lichtenstadt 26) , T. Liska 23) , I. Ludwig 10 ), A. Maggiora 29) , 
M. Maggiora 29 ), A. Magnon 25 ), G.K. Mallot 11 ), IV. Manuilov 24 ), C. Marchand 25 ) , 
J. Marroncle 25 ), A. Martin 28 ), J. Marzec 31) , T. Matsuda 17) , A.N. Maximov 8 ), 
K.S. Medved 8) , W. Meyer 2 ), A. Mielech 28 - 30 ) , Yu.V. Mikhailov 24 ) , M.A. Moinester 26 ) , 
O. Nahle 3) , J. Nassalski 30 ) , S. Neliba 23) , D.P Neyret 25) , V.I. Nikolaenko 24) , 
A. A. Nozdrin 8 ), V.F. Obraztsov 24 ) , A.G. Olshevsky 8 ), M. Ostrick 4 ), A. Padee 31) , 
P. Pagano 28 ), S. Panebianco 25 ) , D. Panzieri 29 - b ) , S. Paul 20 ), H.D. Pereira 10 - 25 ), 
D.V. Peshekhonov 8 ) , V.D. Peshekhonov 8 ) , G. Piragino 29) , S. Platchkov 25 ) , K. Platzer 19 ), 
J. Pochodzalla 16 ), V.A. Polyakov 24 ), A. A. Popov 8 ), J. Pretz 4 ), C. Quintans 15 ), 
S. Ramos 15 ' a ), PC. Rebourgeard 25 ), G. Reicherz 2 ), J. Reymann 10 ), K. Rith 9 ' 11 ), 
A.M. Rozhdestvensky 8 ), E. Rondio 30 ), A.B. Sadovski 8 ), E. Sailer 8 ), V.D. Samoylenko 24 ) , 
A. Sandacz 30 ), M. Sans 19) , M.G. Sapozhnikov 8 ) , LA. Savin 8 ), P. Schiavon 28 ), C. Schill 10 ), 
T. Schmidt 10 ), H. Schmitt 10 ), L. Schmitt 20 ), O.Yu. Shevchenko 8 ) , A. A. Shishkin 8 ), 
H.-W. Siebert 12 - 16 ), L. Sinha 7 ), A.N. Sissakian 8 ), A. Skachkova 29 ) , M. Slunecka 8 ), 
G.I. Smirnov 8 ), F. Sozzi 28 ), V.P Sugonyaev 24 ) , A. Srnka 5 ), F. Stinzing 9 ), M. Stolarski 30 ) , 
M. Sulc 14 ), R. Sulej 31 ), N. Takabayashi 21 ) , V.V. Tchalishev 8 ) , F. Tessarotto 28 ) , 
A. Teufel 9 ), D. Thers 25 ), L.G. Tkatchev 8 ), T. Toeda 21 ), V.I. Tretyak 8 ), S. Trousov 8 ), 
M. Varanda 15 ), M. Virius 23 ), N.V. Vlassov 8 ), M. Wagner 9 ), R. Webb 9 ), E. Weise 3 ), 
Q. Weitzel 20 ), U. Wiedner 19 ), M. Wiesmann 20 ), R. Windmolders 4 ) , S. Wirth 9 ), 
W. Wislicki 30 ), A.M. Zanetti 28 ), K. Zaremba 31 ), J. Zhao 16 ), R. Ziegler 3 ), and 

A. Zvyagin 19 ) 



Universitat Bielefeld, Fakultat fiir Physik, 33501 Bielefeld, Germany f ) 
Universitat Bochum, Institut fiir Experimentalphysik, 44780 Bochum, Germany'-* 
3 ) Universitat Bonn, Helmholtz-Institut fiir Strahlcn- und Kernphysik, 53115 Bonn, Germany') 
Universitat Bonn, Physikalisches Institut, 53115 Bonn, Germany') 
Institute of Scientific Instruments, AS CR, 61264 Brno, Czech Republic 8 ) 

6 ) Burdwan University, Burdwan 713104, India 1 ) 

7 ) Matrivani Institute of Experimental Research & Education, Calcutta- 700 030, India-)' 
Joint Institute for Nuclear Research, 141980 Dubna, Moscow region, Russia 
Universitat Erlangen-Nurnberg, Physikalisches Institut, 91054 Erlangen, Germany f ) 
Universitat Freiburg, Physikalisches Institut, 79104 Freiburg, Germany f ) 
CERN, 1211 Geneva 23, Switzerland 

Universitat Heidelberg, Physikalisches Institut, 69120 Heidelberg, Germany') 

Helsinki University of Technology, Low Temperature Laboratory, 02015 HUT, Finland and University 
of Helsinki, Helsinki Institute of Physics, 00014 Helsinki, Finland 
Technical University in Liberec, 46117 Liberec, Czech Republic 8 ) 

15 ) LIP, 1000-149 Lisbon, Portugal 11 ) 

16 ) Universitat Mainz, Institut fiir Kernphysik, 55099 Mainz, Germany f ) 

17 ) University of Miyazaki, Miyazaki 889-2192, Japan k ) 
18 ' Lebedev Physical Institute, 119991 Moscow, Russia 

19 ) Ludwig-Maximilians-Universitat Munchen, Department fiir Physik, 80799 Munich, Germany f ) 
2 °) Technische Universitat Munchen, Physik Department, 85748 Garching, Germany^ 
21 ) Nagoya University, 464 Nagoya, Japan k ) 

Charles University, Faculty of Mathematics and Physics, 18000 Prague, Czech Republic 8 ) 
23 ) Czech Technical University in Prague, 16636 Prague, Czech Republic 8 ) 

1 : ) State Research Center of the Russian Federation, Institute for High Energy Physics, 142281 Protvino, 
Russia 

CEA DAPNIA/SPhN Saclay, 91191 Gif-sur-Yvette, France 

26 ) Tel Aviv University, School of Physics and Astronomy, 69978 Tel Aviv, Israel 1 ) 

27) ICTP-INFN MLab Laboratory, 34014 Trieste, Italy 
INFN Trieste and University of Trieste, Department of Physics, 34127 Trieste, Italy 
INFN Turin and University of Turin, Physics Department, 10125 Turin, Italy 
Soltan Institute for Nuclear Studies and Warsaw University, 00-681 Warsaw, Poland 111 ) 

31 ) Warsaw University of Technology, Institute of Radioelectronics, 00-665 Warsaw, Poland") 
Also at 1ST, Universidade Tecnica de Lisboa, Lisbon, Portugal 
Also at University of East Piedmont, 15100 Alessandria, Italy 
On leave of absence from JINR Dubna 
Also at Chubu University, Kasugai, Aichi, 487-8501 Japan 
Also at Yamagata University, Yamagata, 992-8510 Japan 
Supported by the German Bundesministerium fiir Bildung und Forschung 
Suppported by Czech Republic MEYS grants ME492 and LA242 

Supported by the Portuguese FCT - Fundacao para a Ciencia e Tecnologia grants 
POCTI/FNU/49501/2002 and POCTI/FNU/50192/2003 
Supported by DST-FIST II grants, Govt, of India 
Supported by the Shailabala Biswas Education Trust 
k ) Supported by the Ministry of Education, Culture, Sports, Science and Technology, Japan; Daikou 
Foundation and Yamada Foundation 

Supported by the Israel Science Foundation, founded by the Israel Academy of Sciences and Human- 
ities 

Supported by KBN grant nr 621/E-78/SPUB-M/CERN/P-03/DZ 298 2000, nr 621/E- 
78/SPB /CERN/P-03 /DWM 576/2003-2006, and by MNII reasearch funds for 2005-2007 
Supported by KBN grant nr 134/E-365/SPUB-M/CERN/P-03/DZ299/2000 



25 



In the nucleoli structure investigations by high energy lepton probes, the region of 
low x corresponds to high parton densities, where new dynamical mechanisms may be 
revealed. The longitudinal structure function gi{x,Q 2 ) is presently the only observable 
which permits the study of low x processes in spin dependent interactions. The existing 
data have been obtained exclusively from fixed-target experiments where the low values 
of x strongly correlate with low values of Q 2 . Therefore theoretical interpretations of the 
results require a suitable extrapolation of the parton ansatz to the low-Q 2 region and 
possibly also an inclusion of nonperturbative mechanisms, which vanish at higher Q 2 . 

Contrary to the spin-independent structure functions, the small-x behaviour of both 
the singlet and the non-singlet part of g\ is controlled by double logarithmic terms, i.e. 
by those terms which correspond to powers of ln 2 (l/x) at each order of the perturbative 
expansion [1]. The double logarithmic effects go beyond the DGLAP evolution and can be 
accommodated in it using special techniques [2-4] . Different approaches permit a smooth 
extrapolation of the obtained g\ to the low-Q 2 region [4, 5] where it may also be comple- 
mented by a non-perturbative component [6]. The double logarithmic terms generate the 
leading small-x behaviour of g± where the relevant Regge poles are expected to have a low 
intercept. 

The region of low x and fixed Q 2 is the Regge limit of the (deep) inelastic scattering 
where the Regge pole exchange model should be applicable. In this model the shape of 
gi at x — > {i.e. at Q 2 <C W 2 where W 2 is the 7*iV centre-of-mass energy squared) is 
parametrised as 

gl(x,Q 2 )~[3(Q 2 ) x -^. (1) 

Here the index i refers to the flavour singlet (s) and nonsinglet (ns) combinations of proton 
and neutron structure functions and aj(0) denotes the Regge trajectory function at zero 
momentum transfer. It is expected that a Sjns (0) < and that a s (0) ~ a ns (0) [7]. This 
behaviour of g\ should translate to the W 2a dependence of the Compton cross-section at 
Q 2 — > where g\ should be a finite function of W 2 , free from any kinematical singularities 
or zeros. 

The spin-dependent structure function of the deuteron gf(x, Q 2 ) has been accurately 
measured in the perturbative region, Q 2 > 1 (GeV/c) 2 [8-12]. Due to the relatively low 
incident energy, the deep inelastic scattering events collected in those experiments cover 
only a limited range of x. The behaviour of g\ at x < 0.001 in the large-Q 2 region is 
unknown due to the lack of data from colliders with polarised beams. 

Measurements at low x and low Q 2 put very high demands on event triggering 
and reconstruction and are very scarce: they were performed only by the SMC at CERN 
on proton and deuteron targets [13]. Here we present new results from the COMPASS 
experiment at CERN on the deuteron longitudinal spin asymmetry A\ and the spin- 
dependent structure function gf in the range 0.001 < Q 2 < 1 (GeV/c) 2 in the photon 
virtuality and 4 • 10~ 5 < x < 2.5 • 10~ 2 in the Bjorken scaling variable. This range is 
essentially the same as that covered by the SMC [13], but the present measurements result 
in about tenfold better precision. They complement our recently published measurements 
obtained in the region 0.004 < x < 0.7 and 1 < Q 2 < 100 (GeV/c) 2 [12]. The data were 
collected during the years 2002 and 2003. They cover the kinematic range presented in 
Fig. [TJ We refer the reader to reference [14] for the description of the 160 GeV/c positive 
muon beam, the two-cell 6 LiD polarised target and the COMPASS spectrometer and to 
Ref. [15] for a detailed description of the analysis. 

The COMPASS data acquisition system is triggered by coincidence signals in ho- 
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doscopes, defining the direction of the scattered muon behind the spectrometer magnets 
and/or by a signal in the hadron calorimeters [16]. Triggers due to halo muons are sup- 
pressed by veto counters installed upstream of the target. COMPASS uses three types of 
triggers: i) inclusive ones, based on muon detection only ii) semi-inclusive triggers, based 
on muon detection and presence of energy deposit in the hadron calorimeters and Hi) 
a calorimetric trigger where only information from the hadron calorimeters is used. The 
low-a: and low-Q 2 region is dominated by semi-inclusive triggers. The contribution of the 
inclusive ones is below 5% for x < 0.001 and exceeds 30% only for x > 0.01. Also the 
contribution of the standalone calorimetric trigger is negligible there. In the kinematic 
region considered here events are characterised by small muon scattering angles and their 
kinematics may be distorted by real photon emission. Therefore in the analysis presented 
here the so-called hadron method [8] is used. This means that all events in our sample 
require the presence of the trajectories of a reconstructed beam muon, a scattered muon 
and at least one additional outgoing particle, together defining an interaction point. The 
presence of hadrons in the final state improves the reconstruction of the interaction point 
and reduces the background of events originating from radiative processes and from the 
muons scattered off atomic electrons. It has been checked that the use of the hadron 
method does not bias the inclusive asymmetries [8]. 

The momentum of the incoming muon, centred around 160 GeV/c and measured in 
the beam spectrometer, is required to be between 140 and 180 GeV/c. The reconstructed 
interaction point has to be located inside one of the target cells. In addition, the extrap- 
olated beam muon trajectory is required to cross entirely both target cells in order to 
equalize the flux seen by each of them. The scattered muon is identified by detectors sit- 
uated behind hadron absorbers and its trajectory must be consistent with the hodoscope 
signals used for the event trigger. 

Events are selected by cuts on the four-momentum transfer squared, Q 2 < 1 (GeV/c) 2 , 
the fractional energy of the virtual photon, 0.1 < y < 0.9, and the scaling variable 
x > 4 ■ 10 -5 . The remaining cuts are the same as those used in the Q 2 > 1 (GeV/c) 2 
analysis [12], with additional quality checks on the interaction point, appropriate to the 
present kinematics [15]. According to the hadron method we also require the most ener- 
getic hadron having > 0.1 (zh is a fraction of the virtual photon energy in the laboratory 
frame, carried by a hadron). 

At low values of x the sample is contaminated by events of muon elastic scattering 
off atomic electrons, /i + e~ — > fi + e~ , occurring at x Me = m e i ectTO n/M = 5.45 -10 -4 (M is the 
proton mass) and at very small scattering angles. To remove such events, cuts are imposed 
on a variable, q9, defined as the product of the angle 9 between the virtual photon and the 
hadron candidate and the sign q of the electric charge of the hadron. Depending upon the 
number of hadron candidates outgoing from the interaction point, the event is rejected if 
— 5<q6<2 mrad or —2 < q6 < mrad depending whether it contains one or two hadron 
candidates 0. The distribution of the qO variable and the x spectrum before and after 
the fie scattering rejection are presented in Fig. [2j The background of fie events which 
remains under the elastic peak is estimated to be smaller than 1% of the data sample. As 
the electromagnetic calorimeter (ECAL) could not be fully used in the present analysis, 
the cuts used for fie scattering rejection presented above are applied in the whole range 
of x to reduce the yield of unwanted radiative events. A study using a small subsample of 
events where the ECAL was available shows that around 50% of those unwanted events 

^ A part of that condition, < q8 < 2 mrad, encompasses misidentified muons and beam halo muons. 
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are excluded from the data sample in this way. The remaining background of radiative 
events accounts for less than about 1% of the data sample. 

The resulting sample consists of 280 million events, out of which about 40% were 
obtained in 2002. This is about 200 times more than in Ref. [13]. The acceptance in the 
(x, Q 2 ) plane after all the cuts is shown by the contour superimposed on Fig. [TJ Average 
values of Q 2 in bins of x are presented in Fig. [3j 

During data taking the two target cells are polarised in opposite directions, so that 
the deuteron spins are parallel or antiparallel to the spins of the incoming muons. The 
spins are inverted every 8 hours by a rotation of the target magnetic field. In 2002 and 
2003 the average beam and target deuteron polarisations were about —0.76 and ±0.51, 
respectively. 

Extraction of the cross-section asymmetry Af in the kinematic region where Q 2 
extends down to about 0.001 (GeV/c) 2 demands special care. The common practise of 
neglecting the m 2 / Q 2 terms ( m M is the muon mass) in the expression for the cross- sect ions 
cannot be applied in this region. Therefore we present below the general spin formalism 
where all the m 2 /Q 2 terms are properly taken into account. The only approximation 
applied is neglecting the m^/E 2 terms (E is the incident muon energy) which are of the 
order of 10~ 7 at our kinematics. In all the formulae we consider the exchange of one 
virtual photon only. The interference effects between virtual Z° and photon exchange in 
the inelastic muon scattering have been measured in an unpolarised experiment [17] and 
found negligible in the kinematic range of current fixed target experiments (see also [18]). 

The polarised inelastic muon-deuteron inclusive scattering cross-section a in the 
one-photon exchange approximation can be written as the sum of a spin-independent 
term a and a spin-dependent term Aa and involves the muon helicity = ±1 



a = a - -h^Aa. (2) 



Eight independent structure functions parametrise the cross-section for a spin-1 target; 
this is twice as many as for the spin- 1/2 case. Apart of the spin-independent structure 
functions F\ and F 2 and the spin-dependent structure functions g\ and g 2 , four additional 
structure functions, b\, b 2 , 63, &4 are needed in the spin-1 case [19]. All these functions 
depend on Q 2 and x. Following previous analyses, cf. Refs. [8,12,13] we neglect &i_ 4 since 
they are predicted to be small [19]. Then the expressions for the cross-sections a and 
Aa and thus the cross-section asymmetries An and A±_ become identical to those for a 
spin- 1/2 target. 

The spin-independent cross-section for parity-conserving interactions can be ex- 
pressed in terms of two unpolarized structure functions Fi and F 2 : 



d 2 a Ana z 

a 



dxdQ 2 Q A x 
where 



xy 2 (l-^)F 1 (x,Q 2 )+ (l- y -ff)F 2 (x,Q" 



(3) 



2Mx VQ 2 
7 = —7= = ~ ( 4 ) 



/Q 2 v 

and v is the energy of the exchanged virtual photon. 

When the muon spin and the deuteron spin form an angle if), the cross-section Aa 
can be expressed as [20] 

Aa = cosip Ami + sinif) cos0 Aa±. (5) 
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Here is the azimuthal angle between the scattering plane and the spin plane. The 
cross-sections Ami and Act_i_ refer to the two configurations where the deuteron spin is 
(anti)parallel or orthogonal to the muon spin; Aery is the difference between the cross- 
sections for antiparallel and parallel spin orientations and Aa± = Actt/ cos0, the dif- 
ference between the cross-sections at angles and + n. The corresponding differential 
cross-sections, which can be written in terms of the two structure functions g\ and gi, are 
given by 



Act,, 



d 2 Acx II \§ixu 2 y 



da;dQ 2 



Q 4 



2 4 



2 2 



Or 



i 2 y 
gi — y 92 



(6) 



and 



d 3 A<T T 8a 2 y f j 2 y 2 



Act t = cos 0Act_i_ = = cos0 — 7 

The relevant asymmetries are 



y 



2m 



2{ 1 + W y)9l + 92 



(7) 



Actn „ Act ± . , 

A\\ = - J A ± = —±. (8) 

2(7 ZCT 

The cross-section asymmetry Ay = (ct^ — ct^)/(ct^ + ct^), for antiparallel (||) and 
parallel spins of the incoming muon and the target deuteron can be obtained from 
the numbers of events iVj collected from each cell before and after reversal of the target 
spins: 

Ni = airmail + P B P T fAf), i = 1, 2, 3, 4, (9) 

where aj is the acceptance, <pi the incoming flux, Tij the number of target nucleons, Pg 
and Pt the beam and target polarisations and / the effective target dilution factor. The 
latter includes a corrective factor p = ct^ 7 / 'ct^°* [21] accounting for radiative events on the 
unpolarised deuteron and a correction for the relative polarisation of deuterons bound 
in 6 Li compared to free deuterons. Average values of / in bins of x for the final data 
sample are presented in Fig. HI The extraction of the spin asymmetry was performed as 
in Refs. [12,22]. The four relations of Eq. corresponding to the two cells and two spin 
orientations, lead to a second-order equation in Ai. This method ensures that fluxes and 
acceptances cancel out in the asymmetry calculation on the condition that the ratio of 
the acceptances of the two cells is the same before and after spin reversal, cf. Ref. [18]. 

The longitudinal and transverse virtual-photon deuteron asymmetries, Af and A d , 
are defined via the asymmetry of absorption cross-sections of transversely polarised photon 
as 

Af = - a T 2 )/{2a T ), A d = (a™ + a^)/{2a T ), (10) 

where ctJ is the 7*-deuteron absorption cross-section for a total spin projection J in the 
photon direction, ctJ l results from the interference between transverse and longitudinal 
amplitudes for J = 0, 1 and a T = (oyj + aj + crJ)/3 is the total transverse photoabsorption 
cross-section. The relation between Af, A d and the experimentally measured Ajf , A^ is 

A d = D(Af + v A d ), A d ± = d(A d 2 -£A d l ), (11) 
where D (the so called depolarisation factor), t], d and £ depend on kinematics: 
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y[(l+ 7 V2)(2-^)-2^/Q 2 ] 
2 / 2 (l-2m2/g2 )(1 + 7 2 ) + 2(l + J R)(l- y - 7 2 2/ 2 / 4 ) ' I J 



7 (1 - y - 7 V/4 - y 2 ml/Q 2 ) 
(l + 7 V2)(l-|//2)-!/X/Q 2 ' 



(13) 



t,/2)(H7V2)-!/W 



A (14) 



7 (l-l//2-y 2 m 2 /Q 2 ) 
£ = 112/0 ' 15 

In view of the small value of rj in our kinematic region the expression for Af in 
Eq. fill I) is reduced to Af ~ A$/D and the possible contribution from the neglected 
term is included in the systematic errors [15]. The virtual-photon depolarisation factor 
D depends on the ratio of longitudinal and transverse photoabsorption cross-sections, 
R = a L /a T . In the present analysis an updated parametrisation of R taking into account 
all existing measurements is used [23] together with an extension to very low values of Q 2 , 
cf. Appendix. Average values of D and R in bins of x are shown in Figs. ??, respectively. 

In order to minimize the statistical error of the asymmetry, the kinematic factors /, 
D and the beam polarisation Pb are calculated event-by-event and used to weight events. 
This approach improves the statistical precision by approximately 8% as compared to 
asymmetry evaluation from events numbers. In the weight calculations a parametrisation 
of Pb as a function of the beam momentum is used. For Pt an average value is used 
for the data sample taken between two consecutive target spin reversals q|. The obtained 
asymmetry is corrected for spin-dependent radiative effects according to Ref. [24] but 
retaining only radiative inelastic tails. 

The final values of Af(x, Q 2 ) are listed in Tabled] with the corresponding average 
values of x and Q 2 . They are also shown as a function of x in Fig. [7J These values confirm, 
with a statistical precision increased by more than an order of magnitude, the observation 
made in Ref. [13] that the asymmetry is consistent with zero for x < 0.01. 

The systematic error of Af contains multiplicative contributions resulting from un- 
certainties on polarisations Pb and Pt, on the dilution factor / and on the function R used 
to calculate the depolarisation factor D. Of these, the largest contribution comes from 
D due to a poor knowledge of R. When combined in quadrature, these errors amount 
to 10-30% (Table [2]). However the most important contribution to the systematic error 
is due to possible false asymmetries which could be generated by instabilities in some 
components of the spectrometer. In order to minimize their effect, the values of Af in 
each interval of x have been calculated for 97 subsamples, each of them covering a short 
period of running time and, therefore, ensuring similar detector operating conditions. An 
upper limit of the effect of the time dependent detector instabilities has been evaluated 
by a statistical approach. Dispersions of the values of Af around their means at each 
value of x were compared with their expected values. Using the Monte Carlo technique 
for a statistical limit estimate [25], values for the false asymmetries were calculated and 
everywhere found to be smaller than the statistical precision. This estimate accounts for 
the time variation effects of the spectrometer components. 



As Pt varies with time, using it in the weight would bias the A\ asymmetry. 



5 



Several other searches for false asymmetries were performed. Data from the two 
target cells were combined in different ways in order to eliminate the spin-dependent 
asymmetry. Data obtained with opposite signs of cell polarisations were compared as 
they may reveal acceptance effects. These searches did not show any significant false 
asymmetry. 

In Fig. [8] results of the present analysis as a function of x are presented together 
with previous measurements by the SMC at 0.01 < Q 2 < 100 (GeV/c) 2 [8, 13]. The 
improvement in the statistical precision at low x is striking. Other data, mostly from the 
deep inelastic scattering region by COMPASS [12], HERMES [11], SLAC E143 [9] and 
SLAC E155 [10], are also presented in Fig. [8] The values of Af, even if originating from 
experiments at different energies, tend to coincide due to the very small Q 2 dependence 
of Af at fixed x. 

The spin dependent structure functions are connected to the virtual photon asym- 
metries in the following way 

^=(il^)W+^)- * = (T&jH + H- (16) 

These formulae are exact; possible contributions from the structure functions &i_4 cancel 
out. Neglecting A\ and making the usual replacement (l + 7 2 )F 2 /(2xFi) = 1 + R, as in the 
spin- 1/2 case and valid if &i_4 = 0, the longitudinal spin structure function gf is obtained 

as 

gf = — — s A t ( 17 ) 

yi 2 x (1 + R) 1 K ' 

The values of gf are listed in the last column of Table [T] and shown in Fig. They 
have been obtained with the F^ parametrisation of Refs. [8,26], cf. Appendix, and with 
the parametrisation of R used in the depolarisation factor. The systematic errors on 
gf are obtained in the same way as for Af, with an additional contribution from the 
uncertainty on F$. Moreover the error of the depolarisation factor was modified. Instead 
of 5D, the error of the quantity D(l + R), 8[D(1 + R)] was considered. The values of 
xgf{x) obtained in this analysis and, for comparison, the SMC [13] and HERMES [11] 
results at Q 2 < 1 (GeV/c) 2 are shown in Fig. [TOl 

The low x data in the kinematic region where W 2 is high and W 2 3> Q 2 , should 
in principle allow testing the Regge behaviour of g\ through its x dependence. These 
conditions are fulfilled by our measurements and thus a fit of Eq. to the gi data from 
the Q 2 range of 0.0025 — 0.25 (GeV/c) 2 in six subintervals of Q 2 ~ const was performed. 
The results of the fit were inconclusive. No information on the singlet intercept, a s (0), 
could be extracted. Thus our data do not provide a test of the Regge behaviour of g\ 
without additional assumptions about its Q 2 dependence. This is due to a limited x 
interval for any given value of Q 2 combined with small measured values of gf. However, 
these data can be compared with models predicting both the x and Q 2 dependence of g\ 
at low values of x and Q 2 [4,6]. A relevant phenomenological analysis is in progress. 

In summary, we have measured the deuteron spin asymmetry Af and its longitudinal 
spin-dependent structure function gf for Q 2 < 1 (GeV/c) 2 over the range 4 • 10~ 5 < 
x < 2.5 ■ 10~ 2 and with a statistical precision more than tenfold better than previous 
experiments. The Af and gf values are compatible with zero for x < 0.01. 
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Appendix 

Knowledge of F 2 (x, Q 2 ) and R(x, Q 2 ) is needed in computations of the dilution fac- 
tor, the radiative corrections, the depolarisation factor and the spin dependent structure 
function gi(x, Q 2 ). It is not sufficient to know these functions only in the kinematic range 
of the analysis since radiative corrections require their knowledge at x > x meas and all 
values of Q 2 including Q 2 = 0, due to radiative "tails". Asymptotic behaviours of F 2 and 
R in the photoproduction limit, Q 2 — >■ 0, are: F 2 ~ Q 2 and R ~ Q 2 (for fixed, arbitrary 
v). These kinematic constraints eliminate potential kinematical singularities at Q 2 = of 
the hadronic tensor defining the virtual Compton scattering amplitude. 

In the analysis, a new SLAC parametrisation of R, Riggs [23], and F 2 parametrisation 
of Ref. [8] is employed. The former, valid for Q 2 > 0.5 (GeV/c) 2 , is extended to lower 
values of Q 2 , including the R ~ Q 2 behaviour at Q 2 = 0, as: 

R(Q 2 < 0.5 (GeV/c) 2 , x) = i?i99 8 (0.5 (GeV/c) 2 , x) ■ (3(1 - exp (-Q 2 /a)) (18) 

with a = 0.2712 and (3 = 1/(1 - exp(0.5/a)) = 1.1880. At Q 2 = 0.5 (GeV/c) 2 the 
function and its first derivative are continuous in the whole x range of our measurements. 
The error on R, 5R, above Q 2 = 0.5 (GeV/c) 2 was taken from Ref. [23] and below Q 2 = 
0.5 (GeV/c) 2 was set to 0.2. For this value of 5R and for the simplest assumption about R 
at Q 2 < 0.5 (GeV/c) 2 and any x (e.g. R = 0.2) there is an approximate agreement (within 
la) with both the value R = at the photoproduction limit and with measurements at 
higher Q 2 from HERA, where R « 0.4 [27]. 

The F 2 of Ref. [8] is valid for Q 2 > 0.2 (GeV/c) 2 and x > 0.0009. At lower values 
of Q 2 and x we used the model of Ref. [26] valid down to Q 2 = and x = 1 -5 and based 
on a concept of generalised vector meson dominance. Two other F 2 parametrisations, 
albeit for the proton [28,29], were also tried together with Ref. [26], to estimate the F 2 
uncertainty, 5F 2 . The former of these parametrisations is based on the parton saturation 
model with recent modifications including the QCD evolution and the latter is a Regge 
motivated fit to all the world data of F 2 , extended into the large Q 2 in a way compatible 
with QCD expectations. They are valid in a range similar to that of Ref. [26]. The SF 2 
uncertainty was taken as the largest difference between the values of the employed F 2 and 
other parametrisations. 
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x range 


(x) 


(Q 2 ) 

[(GeV/c) 2 ] 


(y) 


Af 


gf 


0.000063- 





00004 


0.000052 


0.0068 


0.44 


0.0008 ± 0.0036 ± 0.0034 


0.06 ±0.27 ±0.26 


0.00004 - 





0001 


0.000081 


0.012 


0.49 


-0.0027 ±0.0027 ±0.0017 


-0.22 ±0.23 ±0.14 


0.0001 





00016 


0.00013 


0.021 


0.53 


0.0015 ± 0.0023 ± 0.0014 





13 


± 


0.21 


± 


0.12 


0.00016 - 





00025 


0.00020 


0.034 


0.56 


-0.0007 ± 0.0022 ± 0.0015 


-0 


06 


± 


0.19 


± 


0.13 


0.00025 - 





0004 


0.00032 


0.054 


0.56 


0.0045 ± 0.0022 ± 0.0017 





36 


± 


0.18 


± 


0.14 


0.0004 





00063 


0.00050 


0.085 


0.56 


-0.0022 ± 0.0023 ± 0.0013 


-0 


16 


± 


0.17 


± 


0.09 


0.00063 - 





001 


0.00079 


0.13 


0.55 


-0.0005 ± 0.0025 ± 0.0015 


-0 


03 


± 


0.16 


± 


0.09 


0.001 





0016 


0.0013 


0.20 


0.54 


-0.0035 ± 0.0029 ± 0.0022 


-0 


11 


± 


0.09 


± 


0.09 


0.0016 





0025 


0.0020 


0.32 


0.54 


-0.0023 ± 0.0035 ± 0.0025 


-0 


07 


± 


0.10 


± 


0.07 


0.0025 





004 


0.0031 


0.50 


0.53 


-0.0013 ± 0.0043 ± 0.0034 


-0 


03 


± 


0.10 


± 


0.08 


0.004 





0063 


0.0049 


0.63 


0.43 


-0.0069 ± 0.0061 ± 0.0033 


-0 


11 


± 


0.10 


± 


0.06 


0.0063 





01 


0.0077 


0.68 


0.30 


-0.016 ± 0.010 ± 0.008 


-0 


17 


± 


0.11 


± 


0.09 


0.01 





0158 


0.012 


0.74 


0.20 


0.013 ± 0.019 ± 0.012 





09 


± 


0.13 


± 


0.09 


0.0158 





025 


0.019 


0.82 


0.14 


0.019 ± 0.040 ± 0.019 





09 


± 


0.20 


± 


0.09 



Table 1 : Values of Af and gf with their statistical and systematic errors as a function of x with 
the corresponding average values of x, Q 2 and y. The maximum Q 2 cut is 1 (GeV/c) 2 . Bins in 
x are of equal width in logiox. 



Multiplicative 

variables 

error 


Beam polarisation 


SPb/Pb 


4% 


Target polarisation 


5P t /P t 


5% 


Depolarisation factor 


5D(R)/D(R) 


4 - 30 % 


Dilution factor 


Sf/f 


7% 


Additive 
variables 
error 


Transverse asymmetry 


V ■ SA 2 


< 0.1 • 5A\ tat 


Radiative corrections 


5A? C 


< 0.03 • 5Af at 


False asymmetry 


Afalse 


< 5A\ tat 



Table 2: Decomposition of the systematic error of A\ into multiplicative and additive variables 
contributions. 
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Figure 2: Removal of the /i + e~ — > fi + e~ scattering events. Left: distribution of the variable q9 
(see text for the definition) for events with one (positive or negative) hadron candidate outgoing 
from the primary interaction point. Events between vertical lines are removed from further 
analysis. Note the logarithmic scale on the vertical axis. Right: x distribution of events with 
one negative hadron candidate, before and after \ie event rejection. 
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Figure 3: (Q 2 ) as a function of x for the final data Figure 4: Mean effective dilution factor, (/), as a 
sample. function of x for the final data sample. 
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Figure 5: Mean depolarisation factor, (D), as a Figure 6: Mean values of the ratio R = a L /a T as 
function of x for the final data sample. a function of x for the final data sample. 
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Figure 7: The asymmetry Af(x) as a function of x at the measured values of Q 2 obtained in 
this analysis. Errors are statistical; the shaded band indicates the size of the systematic ones. 
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Figure 8: The asymmetry Af(x) as a function of x at the measured values of Q 2 : the results for 
Q 2 < 1 (GeV/c) 2 obtained in this analysis are compared with previous results at different values 
of Q 2 from COMPASS [12], SMC [8,13], HERMES [11], SLAC E143 [9] and SLAC E155 [10]. 
The E155 data correponding to the same x have been averaged over Q 2 . Errors are statistical. 




Figure 9: The spin dependent structure faction gf(x) as a function of x at the measured values 
of Q 2 obtained in this analysis. Errors are statistical; the shaded band indicates the size of the 
systematic ones. 
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Figure 10: Same as in Fig. [8] but for the quantity xgf. Only data for Q 2 < 1 (GeV/c) 2 are 
shown. 
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